Scaling for QLM communications faster than shannon rate

ABSTRACT

This invention introduces a scaling correction £ to the Quadrature Layered Modulation (QLM) scaling n p  of the communications metrics E b /N o  and SNR. QLM layers n p  communications links over the same frequency and/or path assignment and scales the metrics to the QLM values £n p E b /N o  and £n p   2 SNR in order to maintain the same bit error rate (BER) for all n p . For n p =1 there is no scaling correction £=1. For n p &gt;2 the n p  scaling must be supplemented by a scaling correction £. Identifying £ as a separate parameter whose behavior can be characterized enables one to improve the design and implementation of QLM communications. This introduction of £ enables the bound on QLM data rate performance to become a nearly achievable limit on performance in that one can come arbitrarily close to this bound but can never achieve this bound with each of the n p  QLM layers obeying the Shannon bound.

BACKGROUND OF THE INVENTION

This patent application is a continuation-in-part (CIP) of U.S. Pat. No. 8,917,786.

I. Field of the Invention

The present invention relates to cellular communications and also relates to the Nyquist rate for data symbol transmission, the Shannon bound on communications capacity, Quadrature Layered modulation (QLM), and symbol modulation and demodulation for high-data-rate satellite, airborne, wired, wireless, and optical communications and includes all of the communications symbol modulations and the future modulations for single links and multiple access links which include electrical and optical, wired, mobile, point-to-point, point-to-multipoint, multipoint-to-multipoint, cellular, multiple-input multiple-output MIMO, terrestrial networks, and satellite communication networks. In particular it relates to WiFi, WiFi 802.11ac, WiMax, long-term evolution LTE, 3G, 4G, 5G for cellular communications and satellite communications. WiFi, WiMax use orthogonal frequency division multiplexing OFDM on both links and LTE uses single carrier OFDM (SC-OFDM) on the uplink from user to base station and OFDM on the downlink form base station to user. WiMax occupies a larger frequency band than WiFi and both use OFDM waveforms. SC-OFDM LTE is a single carrier orthogonal waveform version of OFDM which uses orthogonal frequency subbands of varying widths.

II. Description of the Related Art

Bounds on communications capacity are the communications Nyquist rate, the Shannon rate, the Shannon capacity theorem, and the QLM bound. The Nyquist complex sample rate is 1/T=B where B is the signal bandwidth, the Shannon rate W≧B in Hz is the frequency band W=(1+α)B where α is the excess bandwidth required to capture the spillover of the signal spectrum beyond B with a representative value being α=0.25, and the Shannon capacity theorem specifies the maximum data rate C in Bps (bits/second) which can be supported by the communications link signal-to-noise power ratio SNR=S/N over W where S is the signal power and N is the noises power.

The Nyquist rate 1/T is the complex digital sampling rate 1/T=B that is sufficient to include all of the information within a frequency band B over a communications link. Faster than Nyquist rate communications (FTN) transmits data symbols at rates 1/T_(s)≧1/T wherein 1/T_(s) is the data symbol transmission rate in the frequency band B which means T_(s) is the spacing between the data symbols. FTN applications assume the communications links with a data symbol rate equal to the Nyquist rate 1/T_(s)=1/T=B operate as orthogonal signaling with no intersymbol interference (ISI) between the demodulated data symbols. It is common knowledge that the communications data symbol rate 1/T_(s) for orthogonal signaling can be increased to as high as 25% above the Nyquist complex sample rate 1/T_(s)=1/T=B with very little loss in E_(b)/N_(o) and with no perceptible loss in some cases using simple data symbol modulations. Above 25% it has been observed that there is a rapid loss in signal strength. Only a few studies have addressed data symbol rates above 25%. An example in U.S. Pat. No. 8,364,704 transmits digital bit streams at FTN rates and depends on the transmit signal alphabets to be in distinct locations on receive to enable alphabet detection. There is no comparison or proof that the data rate performance is comparable to the performance using conventional orthogonal signaling.

The Shannon bound on the maximum data rate C is complemented by the Shannon coding theorem, and can be defined in equations (1).

Shannon bound and coding theorem

1 Shannon capacity theorem C/W≦log₂(1+SNR)   (1)

2 Shannon coding theorem for the information bit rate R_(b)

-   -   For R_(b)<C there exists codes which support reliable         communications     -   For R_(b)>C there are no codes which support reliable         communications         wherein C in Bps is the channel capacity for an additive white         Gaussian noise AWGN channel in W, “log₂” is the logarithm to the         base 2, and C for the equality “=” is the maximum rate at which         information can be reliably transmitted over a noisy channel         where SNR=S/N is the signal power S to noise power N ratio in W.         It is more convenient to use the communications efficiency C/W         in units of Bps/Hz=b/s/Hz where “b” is the number of information         bits in each communications data symbol. The C/W is the data         rate per unit of bandwith and for convenience of presentation         the C/W will be referenced as the scaled data rate or simply the         data rate. The inequality “≦” reads the C/W cannot exceed the         expression “log₂(1+SNR)”. C/W is almost equal to the expression         with ideal communications links and turbo coding. Ideal         communications links have flat power spectrums over W and are         orthogonal at the Nyquist rate.

MIMO communications enable higher capacities to be supported with multiple independent links over the same bandwidth. This multiple-input multiple-output MIMO requires the physical existence of un-correlated multiple communications paths between a transmitter and a receiver. MIMO uses these multiple paths for independent transmissions when the transmission matrix specifying these paths has a rank and determinant sufficiently large to support the paths being used. In MIMO patent U.S. Pat. No. 7,680,211 a method is disclosed for constructing architectures for multiple input transmit and multiple output receive MIMO systems with generalized orthogonal space-time codes (C₀) which are generalization of space-time codes C and generalizations (H₀) of the transmission matrix (H) that enable the MIMO equation Y=Hf(C,X)+No to be written Y=H₀C₀X+N_(o) which factors out the input signal symbol vector X and allows a direct maximum-likelihood (ML) calculation of the estimate {circumflex over (X)} of X, and wherein Y is the received (Rx) symbol vector, N_(o) is the Rx noise vector, and f(C,X) is a non-separable encoding C of X.

OFDM (orthogonal frequency division modulation) waveform implements the inverse FFT (IFFT=FFT⁻¹) to generate OFDM (or equivalently OFDMA which is orthogonal frequency division multiple access to emphasize the multiple access applications). OFDM uses pulse waveforms in time and relies on the OFDM tone modulation to provide orthogonality. SC-OFDM is a pulse-shaped OFDM that uses shaped waveforms in time to roll-off the spectrum of the waveform between adjacent channels to provide orthogonality, allows the user to occupy subbands of differing widths, and uses a different tone spacing, data packet length, and sub-frame length compared to OFDM for WiFi, WiMax. In addition to these applications the symbol modulations 4PSK (QPSK), 8PSK, 16QAM, 64QAM, 256QAM are used for satellite, terrestrial, optical, and nearly all communication links and with maximum data symbol rates achieved using 256QAM.

SUMMARY OF THE INVENTION

This invention introduces a scaling correction £ for quadrature layered modulation (QLM) which supplements the n_(p) scaling of the signal, and which together are designed to maintain the same error rate performance for all n_(p) being considered. QLM is a FTN layered communications that scales the signal strength to compensate for the inter-symbol (ISI) degradation of the signal strength and to ensure that QLM has the same error rate performance as the original communications link. The parameter n_(p) is the QLM increase in data symbol rate above the current baseband communications data symbol rate for the link being considered, and equivalently is considered to be the number n_(p) of communications layers of QLM. Unless otherwise stated the communications links being addressed for QLM transmit a data symbol rate at the Nyquist rate 1/T_(s)=1/T=B and which means the n_(p) is referenced to the Nyquist rate. Layering of communications over the same link is a preferred way to analyze QLM since it enables one to consider each layer as a separate communications link with the overall capacity being the sum of the component capacities of the n_(p) layers and enables one to prove the scaled Shannon bound applies to each layer.

QLM recognizes that the Shannon bound is a bound on the information capacity supported by the performance metrics E_(b)/N_(o) and SNR=S/N and proves the scaled Shannon bound on data capacity applies to each of the n_(p) layered communications channels. The E_(b) is the data pulse energy per information bit b, N_(o) is the noise power density, S is the signal power, N=N_(o)/T_(s) is the noise power, and SNR is the signal-to-noise power ratio and equivalently is the signal energy to noise energy ratio of the data pulse. Inverting the Shannon bound and using the E_(b)/N_(o) and SNR=S/N definitions enables one to use the Shannon bound to specify the performance metrics C/W, E_(b)/N_(o), SNR as functions of the information bits b supported by each of the data pulses. This is equivalent to specifying that each QLM layer of communications obeys the Shannon bound modified to include the required scaling by n_(p) and £ to maintain the same error rate performance. Communications data rate C/W in units of Bps/Hz=b/s/Hz is the communications efficiency equal to the communications data rate C in Bps per unit of the frequency band W in Hz.

Scaling by n_(p) requires the E_(b)/N_(o) and SNR for the baseband communications link to be increased (scaled) to the values n_(p)E_(b)/N_(o) and n_(p) ²SNR for QLM. For n_(p)=2 this is the exact scaling value and for larger values n_(p)>2 this is a lower bound on the scaling parameter and requires the addition of a “loss” in dB provided by £ to be added to this scaling bound in dB to yield the required QLM scaling loss in dB. This means the required scaled values for QLM are £n_(p)E_(b)/N_(o) and £n_(p) ²SNR.

QLM requires transmitting each layer with a discriminating parameter which enables separation and decoding of each layer. Considering QLM to be an increase in the data symbol rate, the discriminating parameter is the increase in the data symbol rate above the orthogonal Nyquist rate. This increase is detected as a correlated interference by the data symbol detection in the receiver, and enables a demodulation algorithm to unscramble these received and detected correlated data symbol estimates to recover the data symbol information. Discriminating parameters include separate time, frequency, and beam offsets and combinations thereof. A common requirement is that the discriminating parameters provide a means to support a unique inverse so the demodulator can unscramble the detected data symbol estimates

The QLM properties and performance prove that QLM can support data rates at least double the rates supported by the Shannon bound. A new derivation of the bound on the QLM communications capacity is presented which uses a geometric derivation of the Shannon bound modified to apply to QLM. Results are the same as previously derived in QLM patent U.S. Pat. No. 7,391,819 except for the inclusion of £. This geometric derivation provides a new insight into the QLM bound.

Two sets of QLM demodulation algorithms disclosed are the maximum likelihood (ML) data packet and the trellis algorithms. Trellis algorithms are the ML and the co-state maximum a-posteriori (MAP) algorithms. ML algorithms have their roots in the current set of ML convolutional decoding and turbo-decoding algorithms. The MAP algorithms were introduced to provide the correct mathematical technology for these applications and with a decoding and trellis demodulation performance which is at least as-good-as the ML algorithms, These algorithms provide demodulation architectures and implementations of QLM for OFDM, SC-OFDM LTE, MIMO to support at least doubling of the data rates and with applicability to all communications. Performance estimates are the normalized data rate C/W vs. E_(b)/N_(o) and SNR=S/N for the QLM and Shannon bounds, for the reference 4PSK, 8PSK, 16QAM, 64QAM, 256QAM, 4096QAM data symbol modulations, and for representative QLM performance using these symbol modulations. PSK is phase shift keying data symbol modulation and QAM is quadrature amplitude data symbol modulation. Measured bit-error-rate (BER) simulation data is disclosed for trellis demodulation of 4PSK with n_(p)=1-6 layers and for ML data packet demodulation of 4PSK and 256QAM for n_(p)=1-8 layers.

Representative QLM OFDM and QLM SC-OFDM architectures are disclosed for implementation of QLM transmitters and receivers using frequency offset as the differentiating parameter for QLM OFDM and using time offset as the differentiating parameter for QLM SC-OFDM, and implementing the QLM demodulation with ML trellis or MAP trellis algorithms. Frequency layering for QLM OFDM transmits the QLM layers as frequency offset FFT⁻¹ waveforms which are layered over the reference inverse FFT waveform FFT⁻¹ which implements OFDM. For SC-OFDM for LTE the time layering transmits the QLM layers as SC-OFDM waveforms which are offset in time compared to the original SC-OFDM waveform.

A representative MIMO cellular communications link architecture using QLM with MAP trellis demodulation for the transmit and receive signal processing algorithms and supporting block diagrams are developed to illustrate the architecture and implementation,

BRIEF DESCRIPTION OF THE DRAWINGS AND THE PERFORMANCE DATA

The above-mentioned and other features, objects, design algorithms, and performance advantages of the present invention will become more apparent from the detailed description set forth below when taken in conjunction with the drawings wherein like reference characters and numerals denote like elements, and in which:

FIG. 1 describes how to increase the data rate using a pulse waveform.

FIG. 2 describes how QLM increases the data rate for a pulse waveform at a constant frequency bandwidth.

FIG. 3 presents the QLM Lemma 1 which defines the scaling increase of the E_(b)/N_(o) and SNR.

FIG. 4 presents the QLM Lemma 2 bound on QLM C/W for communications and presents the corresponding bounds on the QLM communication link performance metrics C/W, E_(b)/N_(o), SNR.

FIG. 5 defines QLM and lists the QLM properties.

FIG. 6 calculates the time pulse correlation, candidate waveform ψ time response, and the ψ correlation in time.

FIG. 7 calculates the frequency pulse correlation, Discrete Fourier Transform (DFT) waveform frequency response, and the DFT correlation in frequency.

FIG. 8 describes the QLM Maximum Likelihood (ML) symbol packets with n_(s)=2, 3, 4 data symbols in the baseband layer, and illustrates the overlaying of the n_(p)−1 QLM layers on these baseband layers.

FIG. 9 presents the ML data packet measured BER performance vs. the QLM scaled (E_(b)/N_(o))/n_(p) for a pulse waveform with 4PSK modulation for QLM with n_(p)=1,2,4,8 layers for the n_(s)=2,3,4 data symbol packets.

FIG. 10 presents the ML data packet measured BER performance vs. the QLM scaled (E_(b)/N_(o))/n_(p) for a pulse waveform with 256QAM modulation for QLM with n_(p)=1,8 layers for the n_(s)=3,4 data symbol packets, and for a Grey coded 256QAM performance bound.

FIG. 11 describes how the ML and MAP trellis algorithms select the best path j×i from trellis state S_(k−1) to the next trellis state S_(k).

FIG. 12A presents the Matlab code for calculating the best path from trellis state S_(k−1) to trellis state S_(k) and using this selection to update S_(k) for each of the possible mainlobe states.

FIG. 12B presents the Matlab code for calculating the best path from trellis state S_(k−1) to trellis state S_(k) from among the paths stored in S_(k).

FIG. 13A is a flow diagram of a MAP trellis symbol demodulation algorithm.

FIG. 13B is the continuation of the flow diagram of the MAP trellis symbol demodulation algorithm.

FIG. 13C is the continuation of the flow diagram of the MAP trellis symbol demodulation algorithm.

FIG. 14 is a flow diagram of the MAP trellis symbol demodulation algorithm with iteration to include the sidelobes.

FIG. 15 presents the MAP trellis measured BER performance vs. the QLM scaled (E_(b)/N_(o))/n_(p) for a pulse waveform with 4PSK modulation for QLM layers n_(p)=1,2,3,4 using trellis demodulation and for n_(p)=4,6 using a ½-word MAP trellis demodulation.

FIG. 16 calculates the communications data rate C/W versus E_(b)/N_(o) for the new bound on QLM, the Shannon bound for PSK and QAM, and for QLM using ML packet and trellis demodulation with the demodulation parameters in FIG. 18.

FIG. 17 calculates the communications data rate C/W versus SNR for the new bound on QLM, the Shannon bound for PSK and QAM, and for QLM using ML packet and trellis demodulation with the demodulation parameters in FIG. 18.

FIG. 18 lists the MAP trellis and ML demodulation parameters used to calculate the performance expressed as C/W vs. the E_(b)/N_(o) in FIG. 16 and SNR in FIG. 17

FIG. 19 presents the performance data base for C/W vs. E_(b)/N_(o) which demonstrates that the QLM data rate is 2× the Shannon rate for the highest cellular data modulation 256QAM.

FIG. 20 is a representative transmitter implementation block diagram for QLM OFDM.

FIG. 21 is a representative receiver implementation block diagram for QLM OFDM.

FIG. 22 is a representative transmitter implementation block diagram for QLM LTE.

FIG. 23 is a representative receiver implementation block diagram for QLM LTE.

FIG. 24 is a schematic cellular network with the communications link between a base station and one of the users.

FIG. 25 is a representative architecture of this invention disclosure for MIMO systems for application to the communications link in FIG. 24.

DETAILED DESCRIPTION OF THE INVENTION

OFDM, SC-OFDM, LTE applications of quadrature layered modulation QLM in this invention disclosure are illustrated by WiFi and WiMax which uses OFDM on both uplinks and downlinks between the user and base station for cellular communications as well as for communications with satellites, and by LTE which uses OFDM on the downlink and SC-OFDM on the uplink. QLM OFDM adds layers of the OFDM orthogonal data symbol tones at frequency offsets and QLM SC-OFDM adds layers of data symbol tones in separate channels.

FIG. 1 introduces QLM by considering a pulse waveform in the time domain. In 1 the pulse waveform is transmitted at the data symbol rate equal to 1/T_(s)=B where T_(s) is the pulse length, B is the bandwidth, the signal power level 2 is P=A² where “A” is the signal amplitude, and the pulse modulation is phase shift keying PSK with “b” information bits per data symbol. To increase the data symbol rate to n_(p)/T_(s) and the information rate to n_(p)b/T_(s), the pulse waveform is shortened 3 to T_(s)/n_(p) which increases the bandwidth to n_(p)B wherein B=1/T_(s) and requires the transmitted (Tx) power to be increased 4 to P=n_(p)A² in order to keep the same pulse energy per bit 5 is E_(b)=A²T_(s)/b. The corresponding energy-per-bit to noise power ratio 6 is E_(b)/N_(o)=A²/2σ²b where 7 N_(o)=2σ²T_(s) is the noise power density and 2σ² is the “mean square” level of communication noise.

FIG. 2 introduces the QLM scaling law in FIG. 3 QLM Lemma 1 which implements the FIG. 1 increase in the data symbol rate using QLM communications without changing the bandwidth of the pulse waveform by extending the pulses 8 in FIG. 2 over the original pulse length T_(s) and layering these extended data symbol waveforms on top of each other 11 while occupying the same bandwidth b=1/T_(s). The pulse waveforms in each layer 13 have E_(b)/N_(o) values equal to n_(p) times the original E_(b)/N_(o)=A²/2σ²b due to the stretching of each pulse over T_(s) without changing the power level of the pulse. The layers are time synchronized for transmission at ΔT_(s)=T_(s)/n_(p), 2ΔT_(s), . . . , (n_(p)−1)ΔT_(s) offsets 14 respectively for layers 2,3, . . . , (n_(p)−1) relative to the 1^(st) layer at zero offset. This means the signal-to-noise power S/N over B=1/T_(s) is equal to n_(p)^2 times the original S/N due to the addition of the n_(p) pulse power levels 12 over each T_(s) interval and the scaling of E_(b)/N_(o) by n_(p). This scaling of E_(b)/N_(o) in each of the layered communications channels is summarized in equation (2) along with the corresponding scaling of the SNR=S/N over T_(s). We find for the scaling law in FIG. 3 Lemma 1

$\begin{matrix} {\begin{matrix} {{E_{b}/N_{o}} = {{n_{p}\left\lbrack {E_{b}/N_{o}} \right\rbrack}\mspace{14mu}{for}\mspace{14mu}{each}\mspace{14mu}{layer}\mspace{14mu}{or}\mspace{14mu}{channel}}} \\ {= {{n_{p}\left\lbrack {{A^{2}/2}\sigma^{2}} \right\rbrack}/b}} \end{matrix}\begin{matrix} {{S\; N\; R} = {\sum\limits_{n\; p}{b\mspace{11mu}{E_{b}/N_{o}}}}} \\ {= {n_{p}^{2}\left\lbrack {S/N} \right\rbrack}} \\ {= {n_{p}^{2}\left\lbrack {{A^{2}/2}\sigma^{2}} \right\rbrack}} \end{matrix}} & (2) \end{matrix}$ wherein [“o”] is the value of “o” for the baseband communications channel when there is no layering.

FIG. 3 QLM Lemma 1 scaling law derived in (2) for the QLM performance metrics E_(b)/N_(o) and SNR is re-derived from a signal detection viewpoint to complement the derivation of (2) using FIG. 2. QLM received (Rx) demodulation signal processing synchronizes and removes the Rx waveform by performing a convolution of the Rx waveform encoded data symbol with the complex conjugate of this waveform, to detect the correlated data symbols. This convolution is a correlation of the waveform with itself as illustrated in FIG. 6 and FIG. 7 since the waveforms are real and symmetric. Consider a QLM communications link receiver detecting the Tx signal s(t)

$\begin{matrix} {{\cdot {s(t)}} = {\sum\limits_{k}{x_{k}{\psi_{k}\left( {t - {k\;{T_{s}/n_{p}}}} \right)}}}} & (3) \end{matrix}$ where x_(k)ψ_(k) are the data symbol pulses, the x_(k) are the encoded data symbols, and the pulse waveforms ψ_(k) are orthogonal at the Nyquist rate 1/T_(s)=1/T=B, are real and symmetric, are normalized with the square norm equal to unity, and have correlation functions behaving like the triangular correlation function in FIG. 6,7. The Rx signal z(t)=s(t)+v(t) is detected to recover the estimates y_(k) of the Tx data symbols x_(k) by implementing the correlation of the Rx signal with the pulse waveform at the Rx data symbol rate. We find

$\begin{matrix} \begin{matrix} {{y_{k} = {< \psi_{k}}},{{z(t)} >}} \\ {{= {< \psi_{k}}},{{{s(t)} + {v(t)}} >}} \\ {{= {\sum\limits_{n}{< {\psi_{k}\left( {t - {k\;{T_{s}/n_{p}}}} \right)}}}},{{{x_{n}{\psi_{n}\left( {t - {n\;{T_{s}/n_{p}}}} \right)}} + {v(t)}} >}} \\ {= {x_{k} + {\sum\limits_{i}{{c\left( {i} \right)}\left( {x_{k + i} + x_{k - i}} \right)}} + N_{k}}} \end{matrix} & (4) \end{matrix}$ wherein <ψ_(k),z(t)> is the inner product of the waveform ψ_(k) with the Rx signal z(t) after synchronization and translation to baseband, the inner product of ψ_(k) with z(t) is the correlation of the complex conjugate of ψ_(k) with z(t) evaluated at the reference time t=0 for y_(k), correlation coefficients c(|i|) of the waveforms are symmetric c(+i)=c(−i) and equal to the inner product of the waveforms c(i)=<ψ_(k)(t−kT_(s)/n_(p)), ψ_(k+i)(t−(k+i)T_(s)/n_(p))>, and the v(t), N_(k) are the additive white Gaussian noise (AWGN). At the peak the normalized correlation is unity c(0)=1. The first correlation is c(1)=1−1/n_(p), the second is c(2)=1−2/n_(p), . . . , and c(n_(p))=0. The Euclidian distance between the peak correlation value at c(0)=1 and the first correlated signal is equal to 1−c(1)=1/n_(p) and between the second is 1−c(2)=2/n_(p), and so forth. The energy of the first correlated signal is equal to c(1)²(|x_(k+1)|²+|x_(k−1)|²) and behaves as unwanted interference with equal statistical components along and orthogonal to the axis of the detected signal energy E=|x_(k)|²=A²∀k. This means that the squared Euclidian distance between the detected signal x_(k) energy level and the first correlated signal energy level along this detection axis is equal to A²/n_(p) ² for a mainlobe triangular correlation. This requires the signal energy to be increased to n_(p) ²A² to make this shortened distance n_(p) ²(A²/n_(p) ²)=A² equal to the same value A² measured in the absence of any correlation interference for n_(p)=1 and therefore the same error rate performance for n_(p)=2 causing the first correlated interference.. This is a lower bound for n_(p)>2 and is used for the scaling law in (2).

FIG. 5 layered communications property of QLM is demonstrated. We need the QLM FTN architecture to be viewed as a layered communications link to help understand QLM. This layering for integer values of n_(p) is constructed as an overlay of n_(p) independent communications signals s₁(t), s₂(t), . . . , s_(np)(t) in equation (3) and yields

$\begin{matrix} \begin{matrix} {{s(t)} = {\sum\limits_{k}{x_{k}{\psi_{k}\left( {t - {k\;{T_{s}/n_{p}}}} \right)}}}} \\ {= {\sum\limits_{n = 1}^{n\; p}{\sum\limits_{k}{x_{k}{\psi_{k}\left( {t - {\left( {n - 1} \right){T_{s}/n_{p}}} - {k\; T_{s}}} \right)}}}}} \\ {{= {\sum\limits_{n = 1}^{n\; p}{s_{n}(t)}}}{{{where}\mspace{14mu}{s_{n}(t)}} = {\sum\limits_{k}{x_{k}{\psi_{k}\left( {t - {\left( {n - 1} \right){T_{s}/n_{p}}} - {k\; T_{s}}} \right)}}}}} \end{matrix} & (5) \end{matrix}$ with each signal layer transmitting at the orthogonal data symbol Nyquist rate 1/T_(s)=1/T=B. For non-integer values of n_(p) one of the layers has to be transmitted at a different rate. FIG. 2 describes this layering of the communications channels for QLM. QLM is a layered topology for transmitting higher data rates than possible with each layer of communications and is implemented by transmitting each layer with a differentiating parameter which enables separation and decoding of each layer. Each layer or channel has a unique differentiating parameter such as time offset as in (5) or a frequency offset. Each layer or channel obeys the scaled Shannon's laws when using QLM scaling in equations (2).

Shannon bound link metrics for the baseband layer n_(p)=1 are derived by first observing that the Shannon bound in (1) is a bound on the communications performance metrics C/W, Eb/No for n_(p)=1 and 1/T_(s)=1/T=B=W, wherein C/W in Bps/Hz is the capacity per unit of the frequency band W. The Shannon bound specifies the following metrics for 1/T_(s)=1/T=B=W

$\begin{matrix} {{{C/W} = b}\begin{matrix} {{S\; N\; R} = \left( {{s\hat{}b} - 1} \right)} \\ {= {\left( {1/T_{s}} \right)\left( {{A^{2}/2}\sigma^{2}} \right)}} \end{matrix}\begin{matrix} {{E_{b}/N_{0}} = {S\; N\;{R/b}}} \\ {= {\left( {{2\hat{}b} - 1} \right)/b}} \\ {= {\left( {1/T_{s}} \right){\left( {{A^{2}/2}\sigma^{2}} \right)/b}}} \end{matrix}} & (6) \end{matrix}$ which set of metrics is equivalent to the Shannon bound. To prove this we start with the definition 1/T_(s)=1/T=B for n_(p)=1 which proves C/W=(b/T_(s))/(1/T_(s))=b. Inverting the expression for C in (1) and using C/W=b yields SNR=(2^b−1). This enables the definition E_(b)/N_(o)=SNR/b to yield E_(b)/N_(o)=(2^b−1)/b completing the proof. The set of metrics in (6) apply to the baseband communications link for n_(p)=1 which is the first layer of QLM.

FIG. 4 QLM bounds are derived starting with the Shannon bound in (1) expressed as a bound on C/W=log₂(1+SNR). For QLM the SNR is scaled by (2) to read SNR/n_(p)^2 and the scaled Shannon bound applies to a single layer “p” to read (C/W)_(p)=log2(1+SNR/n_(p)^2). The C/W is the sum of the (C/W)_(p) for each of the n_(p) layers and the maximum “max” with respect to n_(p) of the sum C/W over the n_(p) layers yields the QLM bound in (7).

$\begin{matrix} {{{C/W} = {\max\limits_{n\; p}\left\{ {n_{p}{\log_{2}\left\lbrack {1 + {S\; N\;{R/{n_{p}\hat{}2}}}} \right\rbrack}} \right\}}},{B\; p\;{s/{Hz}}}} & (7) \end{matrix}$

FIG. 4 QLM scaling law (2) and bound (7) are re-derived using a geometric approach. A college level geometric derivation of the Shannon bound is modified for QLM to derive the scaling law and bound on C/W for a communication links. The derivation starts by assuming we transmit a log₂(M) bit equal-probable real signal M every T_(s) seconds over the Nyquist frequency band 1/T_(s)=1/T=B=W for complex sampling. The Nyquist sampling theorem specifies there are n=2T_(s)W independent real signal samples (T_(s)W complex samples) in T_(s). Each Tx QLM signal energy in each layer is restricted to an n-dimensional hyperspace of radius [(S/n_(p))(T_(s)/n_(p))^n/2] about the origin where S/n_(p) is the signal power available in each layer and T_(s)/n_(p) is the portion of the T_(s) that is available for demodulation of each layer allowed by the overlap of the signal which is equivalent to the scaling in (2). Each Rx QLM signal energy in each layer is restricted to a radius r=[(S/n_(p))(T_(s)/n_(p))+N T_(s)]^n/2 where the noise energy is by definition N T_(s)=N_(o). We assume a noise energy greater than NT will cause a signal detection error in estimating M. This means the number M of allowable signals that can be Rx without error in the n-dimensional signal space with volume r^n=[(S/n_(p))(T_(s)/n_(p))+NT_(s)]^n/2 is bounded by

$\begin{matrix} \begin{matrix} {M \leq {{\left\lbrack {\left( {{\left( {S/n_{p}} \right)\left( {T_{s}/n_{p}} \right)} + {N\; T_{s}}} \right)/\left( {N\; T_{s}} \right)} \right\rbrack\hat{}n}/2}} \\ {= {{\left\lbrack {{S\; N\;{R/n_{p}^{2}}} + 1} \right\rbrack\hat{}n}/2}} \end{matrix} & (8) \end{matrix}$ which allows the channel capacity C_(p) in Bps in each layer p to be derived

$\begin{matrix} \begin{matrix} {C_{p} = {{\log_{2}(M)}/T_{s}}} \\ {= {W\;{\log_{2}\left\lbrack {1 + {S\; N\;{R/n_{p}^{2}}}} \right\rbrack}}} \end{matrix} & (9) \end{matrix}$ which derives the SNR scaling factor n_(p) ² equivalent to deriving the scaling factor n_(p) for E_(b)/N_(o). The QLM channel capacity C is the sum of the capacity C_(p) in each layer which yields

$\begin{matrix} \begin{matrix} {C = {n_{p}C_{p}}} \\ {= {n_{p}W\;{{\log_{2}\left\lbrack {1 + {S\; N\;{R/n_{p}^{2}}}} \right\rbrack}.}}} \end{matrix} & (10) \end{matrix}$ Re-arranging (10) and taking the maximum with respect to n_(p) yields the QLM bound in 1 in (7) for C/W.

FIG. 4 Lemma 2 replaces the Shannon bound in (1) with the QLM capacity bounds on C/W 1 and E_(b)/N_(o) 2, and performance metrics 4 in (11) for 1/T_(s)=1/T=B=W, and omits the new coding theorem since it has not been addressed in this application. The QLM replacement in (11) for the Shannon bound (1) is derived as follows. The maximum for C/W is derived in (7), the

$\begin{matrix} {\;\begin{matrix} {{{New}\mspace{14mu}{capacity}\mspace{14mu}{bounds}\mspace{14mu} 15\mspace{14mu}{in}\mspace{14mu}{{FIG}.\mspace{14mu} 3}\mspace{14mu}{and}\mspace{14mu}{coding}\mspace{14mu}{theorem}}{{1\mspace{14mu}{C/W}} = {\max\limits_{n\; p}\left\{ {n_{p}\;{\log_{2}\left\lbrack {1 + {S\; N\;{R/{n_{p}\hat{}2}}}} \right\rbrack}} \right\}}}} & \; \\ {{2\mspace{14mu}{E_{b}/N_{o}}} = {\min\limits_{n\; p}\left\{ {\left\lbrack {n_{p}^{2}/\left( {C/W} \right)} \right\rbrack\left\lbrack {{2\hat{}\left( {\left( {C/W} \right)/n_{p}} \right)}{—1}} \right\rbrack} \right\}}} & \; \end{matrix}} & (11) \end{matrix}$

3 New coding theorem

-   -   For R_(b)<C there exists codes which support reliable         communications     -   For R_(b)>C there are no codes which support reliable         communications

4. New performance metrics 26 in FIG. 3

C/W = n_(p)b $\begin{matrix} {{S\; N\; R} = {n_{p}^{2}\left( {{2\hat{}b} - 1} \right)}} \\ {= {\left( {n_{p}^{2}/T_{s}} \right)\left( {{A^{2}/2}\sigma^{2}} \right)}} \end{matrix}$ $\begin{matrix} {{E_{b}/N_{o}} = {S\; N\;{R/\left( {C/W} \right)}}} \\ {= {{n_{p}\left( {{2\hat{}b} - 1} \right)}/b}} \\ {= {\left( {n_{p}/T_{s}} \right)\left( {{A^{2}/2}\sigma^{2}} \right)}} \end{matrix}$ minimum E_(b)/N_(o) in 2 is derived from (7) by solving for E_(b)/No=SNR/(C/W)=n_(p) ²/(C/W)][2^((C/W)/n_(p))−1], in 4 the C/W=n_(p) b/T_(s)W=n_(p)b is the data rate for the n_(p) layers per unit bandwidth of W, the SNR is the scaled value (2) in (6), and the E_(b)/N_(o) is the scaled value (2) in (6). These bounds restate the communications bounds in U.S. Pat. No. 7,391,819 in a more useful format. Wavelet waveform from patent U.S. Pat. No. 7,376,688 reduces the excess bandwidth to α≈0 to meet these requirements for the bounds and performance metrics in addition to the pulse waveforms used to evaluate the QLM performance. The new coding theorem in 3 in equations (11) states that C is the upper bound on the information data rate R_(b) in bits/second for which error correcting codes exist to provide reliable communications with an arbitrarily low bit error rate BER wherein C is defined in 1 in equations (11) and upgrades the Shannon coding theorem 2 in equations (1) and new data symbol rate n_(p)/T_(s) is n_(p) times the Nyquist rate 1/T=B for a single channel.

FIG. 5 defines QLM to be a layered communications link equivalent to increasing the data symbol rate with the derived properties for the layering in (5), the Lemma 1 scaling in (2), Lemma 2 bounds in (11), and the Lemma 2 link metrics in (11).

FIG. 6 in time offset units of T_(s) calculates the ideal triangular correlation 10, an example waveform 11 designated by ψ, and the waveform ψ correlation 12. Parameters of interest for this example square-root raised-cosine waveform are the waveform length L=3, M=16, and excess bandwidth α=0.22, the mainlobe 13 which extends over a 2T_(s) interval, and the sidelobes 14 which fall outside of the mainlobe. Parameter L is the waveform length in units of M=16 samples and M is the number of digital samples between adjacent waveforms at a Nyquist symbol rate=1/T_(s) for which T_(s)=MT where 1/T is the digital sample rate and α is a measure of the roll-off of the frequency response. The ideal triangular correlation is the correlation for the pulse waveform of length T_(s) in FIG. 2 and FIG. 7 demonstrates that for waveforms of interest for QLM the triangular correlation approximates the mainlobe correlations for QLM waveforms.

FIG. 7 in frequency offset units of 1/T_(s) presents a N=64 point discrete Fourier transform DFT for OFDM in frequency, the correlation function, and an ideal correlation function in frequency for a pulse waveform in frequency. The correlation function closely approximates the waveform and the ideal correlation with a triangular mainlobe, and has a mainlobe 22 extending over 2/T_(s)=2B Hz interval with low sidelobes 23 wherein the symbol rate 1/T_(s) is at the Nyquist rate and equal to the bandwidth B for each channel

QLM demodulation algorithms include ML and MAP trellis algorithms, ML data packet algorithms, sequential relaxation algorithms, equalization, super-resolution and other algorithms. Sequential relaxation algorithms process the stream of data symbol samples to find the best estimates of the data symbols using sequential relaxation methods. These are computationally intensive and the least developed of the other algorithms. We first consider the ML algorithms in the architecture of data packets.

FIG. 8 presents a representative set of ML QLM data packet architectures for time differentiating parameters which support implementable ML demodulation algorithms with acceptable demodulation losses and with applicability to a broad spectrum of communications systems. ML algorithms solve the QLM system equations for the correlated signals to find estimates of the data symbols in each of the QLM layers and are the least complex of the demodulation algorithms, are limited to a relatively few data symbols, and require the waveform correlation functions to approximate the ideal pulse correlations in FIG. 6,7. The three data packet architectures in FIG. 8 are a n_(s)=2-symbol group 101, a n_(s)=3-symbol group 102, and a n_(s)=4-symbol group 103, with each group respectively consisting of n_(s)=2,3,4 contiguous symbols 104,106,108 for the first QLM layer n_(p)=1. QLM layers for n_(p)>1 are overlayed on each group of symbols such that the additional QLM symbols stay within the composite bandwidth for each symbol group. These architectures assume the sidelobes of the data symbols will be included in the trellis algorithms by an iterative algorithm similar to the sidelobe algorithm disclosed in FIG. 14. The QLM ML packet demodulation algorithm finds an estimate {circumflex over (X)} of the Tx data symbol vector X in a QLM data packet in FIG. 8 by solving the Rx vector equation (12) in the detected data symbol Rx vector Y. Each row “k” of the Rx vector equation (12) is the detection equation (4) which defines the detected signal y_(k) as a linear function of the correlated subset of Tx data symbols {d_(k)|∀k} within the packet and the AWGN. The correlation coefficients {c(|i|)|i=0, . . . , N_(x)−1} in (4) are the elements of the Rx system matrix H for the mainlobe correlation which is a symmetric N_(x) by N_(x) toeplitz correlation matrix over the ML data symbol packet in FIG. 9 where N_(x) is the number of packet symbols N_(x)=[(n_(s)−1)n_(p)+1] and n_(s)=2,3,4 is the number of baseband symbols in a packet. Elements H(k,n) of H are the correlation values of the data symbol waveforms H(k,n)=<ψ(t−kT/n_(p)), ψ(t−nT/n_(p))>=c(|k−n|) with diagonal elements H(k,k)=c(0)=1∀k where indices “k,n” are over the data symbols in the data packet. The size of the packet is limited to n_(s)=4 baseband data symbols to prevent the determinant of H from becoming too small and thereby causing the dispersion of the inverse H⁻¹ of H to be too large for reliable estimation of the data symbols in AWGN. The ML waveforms used to measure the performance data are the square pulses in FIG. 8 with triangular correlations in FIG. 6,7 and no sidelobes.

A ML data symbol packet (pulse group) in FIG. 8 consists of n_(s)=2, 3, 4 baseband square data symbols of length T_(s), 1/T_(s) for time, frequency square pulse waveforms. QLM overlays layers of data symbol pulses over each pulse group n_(s)=2, 3, 4 and the matched filter detection in (2) recovers the elements of the Rx vector Y=[y₁ y₂, . . . , y_(Nx)]^(T) in (12) of the Tx data symbol vector X=[x₁, x₂, . . . , x_(Nx)]^(T) where N is the AWGN noise vector and the elements {x_(k)},{y_(k)} are defined in (4). From (4) the Rx vector equation for the mainlobe correlation is Y=HX+N  (12)

where

-   -   Y=N_(x)×1 detected symbol vector     -   H=correlation matrix of data symbols     -   X=data symbol vector over layered Filters/channels     -   N=demodulation plus link noise vector

$\begin{matrix} {N_{x} = {{number}\mspace{14mu}{of}\mspace{14mu}{QLM}\mspace{14mu}{channels}\mspace{14mu}({filters})}} \\ {{\left( {n_{s} - 1} \right)n_{p}} + 1} \end{matrix}$ wherein the Rx data symbol vector Y has components {y_(i)} over the set of channels or filters indexed on “i”, Tx data symbol vector X has components {x_(i)} indexed over the same set of channels or filters, and Rx noise vector N has components {N_(i)} indexed over the same set of channels or filters. The ML demodulation algorithm is derived in 4 in equations (13) for the system equations (12). We find

$\begin{matrix} {{M\; L\mspace{14mu}{demodulation}\mspace{14mu}{algorithm}}\begin{matrix} {{1\mspace{14mu} M\; L\mspace{14mu}{Cost}\mspace{14mu} J} = {\left\lbrack {Y - {H\; X}} \right\rbrack^{\prime}{Q^{- 1}\left\lbrack {Y - {H\; X}} \right\rbrack}}} \\ {= {( - )\mspace{11mu}{exponent}\mspace{14mu}{of}\mspace{14mu}{the}\mspace{14mu} M\; L\mspace{14mu}{probability}\mspace{14mu}{densisty}}} \\ {function} \end{matrix}} & (13) \end{matrix}$

-   -   2. ML solution {circumflex over (x)} minimizes J         {circumflex over (x)}=[H′ Q ⁻¹ H] ⁻¹ H′ Q ⁻¹ Y     -   3 The noise covariance Q is

$\begin{matrix} {Q = {E\left\lbrack {N\; N^{\prime}} \right\}}} \\ {= {2\sigma^{2}H}} \end{matrix}$ since the noise is correlated by the same filter overlap as the signal. This reduces the ML solution 2 to:the reduced ML demodulation algorithm {circumflex over (x)}=H ⁻¹ Y  (14) wherein H′ is the conjugate transpose of H, the inverse H⁻¹ exists for all waveforms of interest, and “2σ²” is the root-mean-square rms data symbol detection noise.

The ML algorithm for QLM in equations (12),(13),(14) equally applies when using partial symbol ΔT_(s) integration wherein the integration intervals are over the separation between the consecutive layers of communications so that ΔT_(s)=T_(s)/n_(p). The correlation matrix for a full symbol integration is a symmetric Toeplitz matrix whereas using a ΔT_(s)=T_(s)/n_(p) integration interval reduces the matrix to a triangular form thereby offering a decrease in computational complexity and a more robust solution since the determinant of the correlation matrix will not become too small with relatively large data packets. However, the requirement to demodulate with a matched waveform makes it extremely difficult to use partial symbol ΔT_(s) integration.

FIG. 9 measures the ML bit error rate BER performance for uncoded 4PSK pulse waveforms for no QLM wherein n_(p)=1 and for n_(p)=1,2,4,8 layers of of QLM modulation using n_(s)=2,3,4 data symbol groups and implementing the ML symbol demodulation algorithm in (14). Performance validates QLM technology and scaling. QLM performance is plotted as bit error rate BER versus the QLM normalized value (E_(b)/N_(o))/n_(p) of the E_(b)/N_(o). Measured performance values are from a direct error count Monte Carlo simulation of the trellis algorithm and are plotted as discrete measurement points.

FIG. 10 measures the ML BER performance for uncoded 256QAM pulse waveforms for no QLM wherein n_(p)=1 and for n_(p)=8 layers of QLM modulation using n_(s)=3,4 data symbol groups and implementing the ML symbol demodulation algorithm in (14). The 256QAM bound assumes a Grey code assignment of the data bits to the 256 data points in the complex plane and therefore provides improved performance compared to the 256 point regular data mapping used in the simulations.

FIG. 11 introduces the trellis path selection in the ML and MAP trellis demodulation algorithms which unscramble the detected correlated signals in (4) to recover the information in the Tx data symbols. After down-conversion and synchronization the Rx QLM signal z(t) is processed as described in (4) to strip off the waveform ψ_(k) by convolving the z(t) with the complex conjugate of ψ_(k), to recover the correlated data symbol y_(k) estimate of the Tx data symbol x_(k). Trellis processing of the y_(k) recovers the m_(k)=1,2, . . . , 2^b_(s) modulation states (alphabet elements) for each x_(k) wherein b_(s) is the number of uncoded bits in each data symbol x_(k), the information bits b=R b_(s) wherein R is the code rate, and the modulation states m_(k) are the encoded complex symbols representing x_(k).

This convolution is a correlation of the waveform with itself as illustrated in FIG. 6, since the waveforms are real and symmetric. QLM trellis demodulation algorithms process the mainlobes of the waveform correlation algorithms in the first pass and pick up the impacts of the sidelobes in subsequent passes. This iterative architecture provides a means to demodulate the stream of over-lapping and correlated QLM data symbol waveforms with reasonable computational requirements by processing the mainlobes and including the sidelobes in subsequent iterations. Correlation data in FIG. 6, 7 suggest this approach since the sidelobes are relatively low compared to the mainlobes for the bandwith efficient waveforms being addressed.

MAP trellis algorithms have the correct architecture for reducing the number of trellis paths to enable higher order modulations to be used for larger numbers of overlapped QLM channels, compared to the maximum likelihood ML trellis algorithms which are derived from the ML trellis algorithms used for convolutional and Turbo decoding. Another critical property disclosed in patent U.S. Pat. No. 7,337,383 is that the MAP trellis algorithms are at least-as-good-as the ML trellis algorithms.

FIG. 12 trellis paths j×i for the ML and MAP trellis algorithms are paths from the set of xi nodes defining the previous state S_(k−1) to the set of jx nodes defining the current state S_(k). Path j×i reads from left to right, “j” is the new data symbol state, “i” is the last data symbol state, “x” is the string of path data symbol states between “j” and “i”, and “k” is step k corresponding to the detected symbol y_(k) in (4). Nodes designated by circles in FIG. 12 are the sets of symbols identified by their indices, state S_(k−1) is the set of xi nodes for state k−1 and state S_(k) is the set of jx nodes for state k. The data symbols in the transition paths j×i between nodes are the mainlobe correlation values in FIG. 6,7 encoded with the data modulations m_(x). For a QLM signal which has n_(p) simultaneous communications channels, there are (2n_(c)−1) correlation values over the mainlobe which means each state S_(k−1), S_(k) has 2^b_(s)(2n_(p)−2) nodes identified as xi,jx respectively and there are n_(s)^(2n_(p)−1)=2^b_(s)(2n_(c)−1) j×i paths from state S_(k−1) to state S_(k) data symbols wherein 2^b_(s) is the number of possible modulaton symbols m_(x) or un-encoded states of each data symbol, Sidelobe correlations for n_(c)>n_(p) are included in iterative MAP trellis demodulation algorithms and wherein n_(c) is the number of one-sided correlation values including the sidelobes. As illustrated in FIG. 12 the ML algorithm selects the best path 25 by finding the xi which provides the best transitioning metric j×i from xi to jx from among the possible choices 26 for xi. The MAP algorithm selects the best path 27 by finding the jx which provides the best transitioning j×i from xi to jx from among the possible choices 28 for jx. The MAP is a co-state formulation of the trellis algorithm since it replaces the search over jx with a search over xi thereby enabling the number of paths specified by xi to be reduced along with the number of states in the trellis diagram and the computational complexity, and with minimal impact on performance.

MAP and ML trellis algorithms implement the decisioning rules in (15) for updating the path and state metric α_(k)(S_(k)) for each node in the trellis diagram

$\begin{matrix} \begin{matrix} {{\alpha_{k}\left( S_{k} \right)} = {\min\left\{ {{{\alpha_{k - 1}\left( S_{k - 1} \right)} + {R_{k}\left( {S_{k - 1}->S_{k}} \right)}} = {R_{k}({jxi})}} \right\}}} \\ {= {\text{new metric for state}\text{k}}} \end{matrix} & (15) \end{matrix}$ wherein

min   is  with  respect  to  i   for  M L min   is  with  respect  to  j   for  M A P $\begin{matrix} {{\alpha_{k}\left( S_{k} \right)} = {p\left( y_{k} \middle| S_{k} \right)}} & {{for}\mspace{14mu} M\; L} \\ {= {p\left( S_{k} \middle| y_{k} \right)}} & {{for}\mspace{14mu} M\; A\; P} \end{matrix}$ $\begin{matrix} {{R_{k}\left( {S_{k - 1}->S_{k}} \right)} = {R_{k}\left( {j\; x\; i} \right)}} & \; \\ {= {{{{y_{k} - {\hat{y}}_{k}}}^{2}/2}\sigma^{2}}} & {{for}\mspace{14mu} a\mspace{14mu} M\; L\mspace{14mu}{metric}} \\ \left. {= {{{x_{k}}\hat{}2} - {2*{{Re}\left( {y_{k}{\hat{y}}_{k}} \right)}*}}} \right\} & {{for}\mspace{14mu} a\mspace{14mu} M\; A\; P\mspace{14mu}{metric}} \end{matrix}$ wherein (o)* is the complex conjugate of (o), the choice of a “min” or a “max” depends on the definition of R(j×i), the y_(k) refers to the correlated data symbol measurement at state or step k, the estimate ŷ_(k) of y_(k) is the correlated sum of the estimates {{circumflex over (x)}_(k)} of the data symbols {x_(k)}, and the R_(k)(S_(k−1)−>S_(k))=R_(k)(j×i) is the state k transition decisioning metric for the transition written symbolically as S_(k−1)−>S_(k) and realized in Matlab notation as “j×i”.

FIG. 12 lists the Matlab code in a pseudo-code format for the calculation of the best trellis paths in and for the state updates for the ML and MAP trellis algorithms. For the best trellis path, the code translates the trellis algorithms in (15) written as a state metric update equation, into a Matlab code format for calculation of the state S_(k) metric α_(k)(S_(k)) for the best trellis transition path. ML trellis algorithm finds the last symbol “i” with the best path “jx” to the new symbol “j” whereas the MAP algorithm finds the new symbol “j” with the best path “xi” from the last symbol “i” to support a reduction in the number of paths “xi” when necessary to reduce the computational complexity. For the state update for k≧D the state metric S_(k) is upgraded for this new path j×i by the update operations S_(k)(:,jx)=[y_(k)(j×i); S_(k−1)(1:D−1, xi)], S_(k)(:,xi)=[y_(k)(j×i); S_(k−1)(1:D−1, jx)] respectively for ML, MAP using Matlab notation which replaces the column jx, xi vector with the column xi, jx vector after the elements of xi, jz have been moved down by one symbol and the new estimated symbol {circumflex over (x)}_(k) added to the top of the column vector which is the row 1 element.

FIG. 13A is a summary implementation flow diagram of a representative MAP trellis data symbol demodulation algorithm with error correction decoding and with the option of including the sidelobes in subsequent iterations. Basic elements of the algorithm are the trellis parameters 170, MAP trellis algorithm 171 structured by the parameter set in 170, and followed by trellis error correction decoding 172 which also could be combined with the trellis demodulation, and with the option 175 to include the sidelobes in an iteration(s) or to end the iteration(s) and handover 174 the data estimates to the receiver. The parameter set 170 is applicable to the mainlobes of the data symbol waveforms with two examples illustrated in FIG. 6,7. The correlation function for the data symbol waveforms defined in FIG. 6,7 have low sidelobes compared to the mainlobes which suggests it is sufficient to consider mainlobe correlation in the MAP trellis algorithm with the sidelobes included in subsequent iterations. A one-sided correlation vector C is used since the correlation is symmetrical about the origin, whose elements are the mainlobe correlation coefficients C=[c(0), c(1), . . . , c(n_(p)−1)] wherein c(0)=1 with waveform normalization. Waveform modulations emphasized are 4PSK, 8PSK, 16QAM, 64QAM, 256QAM with b_(s)=2, 3, 4, 6, 8 uncoded bits per data symbol and n_(s)=2^b_(s)=4, 8, 16, 64, 256 states per data symbol.

FIG. 13B,13C is a detailed implementation flow diagram of the MAP trellis data symbol demodulation algorithm in FIG. 13A for mainlobe demodulation. FIG. 13B uses the trellis parameters 170 in the initialization 150 prior to the start k=0 with an empty shift register (SR) which normally contains the 2n_(p)−1 correlated data symbols for the mainlobe for each node being addressed. At k=0 the path metric α₀(xi) is set equal to a negative or zero initialization constant for all of the nodes xi=0,1,2, . . . , (n_(s)^(2n_(p)−2)−1} of the trellis diagram where α₀(xi) is the logarithm of the state S₀ path metric at k=0 for node xi, the receive Rx symbols are indexed over k with k=0 indicating the initial value prior to the Rx symbol k=1, nodes of the trellis diagram are the states of the shift register SR, and state S_(k) refers to the trellis diagram paths and metrics at symbol k in the trellis algorithm.

In FIG. 13B loop 151 processes the Rx symbols k=1,2, . . . , n where the index k also refers to the corresponding algorithm steps and the states of the trellis algorithm. In 152 the Rx signals are pulse detected to remove the carrier frequency and waveform to recover the normalized correlated data symbol y_(k). For each Rx symbol y_(k) the state transition decisioning metrics R_(k)(j×i) are calculated by the transition metric equations {R_(k)(j×i)=−|y_(k)−ŷ_(k)(j×i)|^2} for a ML metric, {R_(k)(j×i)=|ŷ_(k)|^2−2Real(y_(k)ŷ_(k)(j×i))*} for a MAP metric wherein (o)* is the complex conjugate of (o), for {R_(k)(j×i)=Metric(y_(k),ŷ_(k))} for another metric, for all possible transition paths {j×i} from the previous state S_(k−1)(xi) at node xi to the new state S_(k)(jx) at node jx in the trellis diagram and where ŷ_(k)(j×i) is the hypothesized normalized detected correlated symbol k for the path j×i. For a mainlobe correlation function the ŷ_(k)(j×i) is defined by the equation ŷ_(k)(j×i)=c(n_(p)−1)[sr(1)+sr(2n_(p)−1)]+ . . . +c(1)[sr(n_(p)−1)+sr(n_(p)+1)]+c(0)[sr(n_(p))] which calculates ŷ_(k)(j×i) as the correlated weighted sum of the elements of the shift register SR=[sr(1),sr(2), . . . , sr(2n_(p)−1)]′ with {circumflex over (x)}_(k)=sr(n_(p)), {circumflex over (x)}_(k−1)=sr(n_(p)−1), {circumflex over (x)}_(k+1)=sr(n_(p)+1), . . . where c(0)=1, the normalized data symbol estimates {{circumflex over (x)}_(k)} correspond to the transition index j×i, and the state k estimated symbol {circumflex over (x)}_(k) is the SR center element sr(n_(p)) with correlation coefficient c(0)=1. Symbols move from left to right starting with “j” with each new Rx symbol or step in the trellis recursion algorithm, and ending with “i”. With this convention “j” is indexed over the states of sr(1), “x” is indexed over the current states of sr(2), . . . , sr(2n_(p)−2), and “i” is indexed over the states of sr(2n_(p)−1). Index over the paths of the trellis diagram is defined by the equation j×i=sr(1)+n_(s) sr(2)+n_(s)^2 sr(3)+ . . . +n_(s)^(2n_(p)−2) sr(2n_(p)−1)−1=0,1,2, . . . , n_(s)^(2n_(p)−1)−1 when the contents of the SR elements are the indices corresponding to the assumed data symbol state values.

In FIG. 13B loop 153 calculates the best trellis transition paths from state S_(k−1) to the new state S_(k) for xi=0,1,2, . . . , n_(s)^(2n_(p)−2)−1. In 154 the path metric α_(k)(S_(k)) is defined by the recursive logarithm equation α_(k)(S_(k))=α_(k−1)(S_(k−1))+R_(k)(S_(k−1)−>S_(k)) which can be rewritten as α_(k)(xi)=α_(k−1)(jx)+R_(k)(j×i) since α_(k)(S_(k))=p(S_(k)|y_(k))=α_(k)(xi), α_(k−1)(S_(k−1))=p(S_(k−1)|y_(k−1))=α_(k−1)(jx), R_(k)(S_(k−1)−>S_(k))=R_(k)(j×i) from equation (5). The best path metric α_(k)(xi) for each new node jx is chosen by the decisioning equation α_(k)(xi)=min{α_(k−1)(jx)+R_(k)(jx)} with respect to the admissible “j”. For each jx, the corresponding xi yielding the highest value of the path metric α_(k)(xi) is used to define the best path and new symbol {circumflex over (x)}_(k−D).

FIG. 14 is a flow diagram of an iterative version of the MAP data symbol trellis algorithm in FIG. 13 wherein the iteration is used to incorporate the effects of the sidelobes of the correlation function C into the trellis state transition metric R_(k)(j×i) in 152 in FIG. 13B with a relatively small increase in computational complexity. In FIG. 14 the algorithm is initialized with the parameter set in 176 which is the parameter set in 170 in FIG. 13A with the partitioning of the correlation vector C into the mainlobe vector C₀ plus the sidelobe vector C₁ and adding the specification of the stopping rule for the iterations. For a correlation vector C=[c(0),c(1), . . . , c(n_(p−)1),c(n_(p)), . . . ,c(n_(c))] consisting of n_(c) correlation coefficients, the mainlobe vector is C₀=[c(0),c(1), . . . , c(n_(p)−1),0,0, . . . , 0] and the sidelobe vector is C₁=[0,0, . . . , 0,c(n_(p)), . . . , c(n_(c))] to partition C=C₀+C₁. The algorithm starts 177 by implementing the MAP algorithm in FIG. 12.

In FIG. 14 the output data symbol estimates 164 in FIG. 13C are the outputs from the MAP algorithm in 177 and are used to calculate the a-priori estimated sidelobe contribution ŷ_(k|1) in 179 to ŷ_(k)(j×i) which is the hypothesized normalized detected correlated symbol k for the path j×i in the calculation of the metric R_(k)(j×i) in 152 in FIG. 13B. In this implementation 180 of the trellis algorithm 171 in FIG. 14A, the ŷ_(k)(j×i)=ŷ_(k|0)(j×i)+ŷ_(k|1) is the sum of the hypothesized mainlobe contribution ŷ_(k|0)(j×i) using the non-zero C₀ coefficients as described in FIG. 13A and the sidelobe contribution ŷ_(k|0) using the data symbol estimates from 177 and the non-zero C₁ coefficients and wherein the subscripts “k|0” reads “index k given C₀” and “k|1” reads “index k given C₁”. From 152 in FIG. 13B we find the ŷ_(k|0)(j×i) is defined by the equation ŷ_(k|0)(j×i)=c(n_(p)−1)[sr(1)+sr(2n_(p)−1)]+ . . . +c(1)[sr(n_(p)−1)+sr(n_(p)+1)]+c(0)[sr(n_(p))] which calculates ŷ_(k|0)(j×i)=ŷ_(k)(j×i) in FIG. 13B as the correlated weighted sum of the elements of the shift register SR=[sr(1),sr(2), . . . , sr(2n_(p)−1)]′ with {circumflex over (x)}_(k)=sr(n_(p)), {umlaut over ({circumflex over (x)}_(k−1)=sr(n_(p)−1), {circumflex over (x)}_(k+1)=sr(n_(p)+1),. . . where c(0)=1, the normalized data symbol estimates {{circumflex over (x)}_(k)} correspond to the transition index j×i, and the state k estimated symbol {circumflex over (x)}_(k) is the SR center element sr(n_(p)) with correlation coefficient c(0)=1. Symbols move from left to right starting with “j” with each new Rx symbol or step in the trellis recursion algorithm, and ending with “i”. With this convention “j” is indexed over the states of sr(1), “x” is indexed over the current states of sr(2), . . . ,sr (2n_(p)−2), and “i” is indexed over the states of sr(2n_(p)−1). Index over the paths of the trellis diagram is defined by the equation j×i=sr(1)+n_(s) sr(2)+n_(s)^2 sr(3) + . . . +n_(s)^(2n_(p)−2) sr(2n_(p)−1)−1=0,1,2, . . . , n_(s)^(2n_(p)−1)−1 when the contents of the SR elements are the indices corresponding to the assumed data symbol state values. The sidelobe contribution is equal to ŷ_(k|1)=c(n_(p))({circumflex over (x)}_(k−n) _(p) +{circumflex over (x)}_(k+n) _(p) )+c(n_(p)+1)({circumflex over (x)}_(k−1−n) _(p) +{circumflex over (x)}_(k+1+n) _(p) )+c(n_(p)+2)({circumflex over (x)}_(k−2−n) _(p) +{circumflex over (x)}_(k+2+n) _(p) )+ . . . until the end of the sidelobe correlation coefficients or the end of the data symbol bit estimates and wherein {circumflex over (x)}_(k−n) _(p) is the data symbol estimate in 180 for symbol k−n_(p). A stopping rule in 182, 184 is used to decide if another iteration is required. When another iteration is required the data symbol bit estimates are used to update the calculation of the a-priori contribution ŷ_(k|1) of the sidelobes to the ŷ_(k)(j×i) in the modified trellis algorithm 184. After the iteration is stopped the outputs are error corrected decoded 183 and handed over to the receiver. A possible alternative path is to use the decoded output estimates 185 to re-start 184 the MAP trellis demodulation algorithm with sidelobes.

FIG. 15 measures the trellis decoding performance for uncoded 4PSK pulse waveforms for no QLM wherein n_(p)=1 and for n_(p)=2,3,4 layers of QLM modulation using the trellis data symbol decoding algorithm in FIG. 13 and for n_(p)=4,6 layers by implementing a ½-word trellis data symbol decoding algorithm. Performance validates QLM technology and scaling. QLM performance is plotted as BER versus the QLM normalized value (E_(b)/N_(o))/n_(p) of the E_(b)/N_(o). Measured performance values are from a direct error count Monte Carlo simulation of the trellis algorithm and are plotted as discrete measurement points. It is expected that higher order modulations 8PSK, 16QAM, 64QAM, 256QAM will have the same relative trellis coding performance for n_(p)=1,2,3,4,5,6 layers of communications.

FIG. 16,17 calculate the QLM data rate C/W vs. E_(b)/N_(o), SNR respectively using the MAP trellis demodulation and ML packet demodulation algorithms. Performance calculated are the Shannon bound, new bound, 4PSK, 8PSK, 16QAM, 64QAM, 256QAM, 4096QAM, MAP trellis demodulation, and ML demodulation communications. Equation (1) is used to calculate the Shannon bound, equation (7) is used to calculate the new bound, modulation performance for 4PSK, 8PSK, 16QAM, 64QAM, 256QAM, 1024QAM assumes turbo coding performance provides a performance almost equal to the Shannon bound and with information rate b=1.5,2,3,4,6,8 bits=Bps/Hz respectively assuming coding rates R=¾, ⅔, ¾, ⅔, ¾, ⅔ and bits per symbol b_(s)=1.5,2,3,4,6,8 bits respectively for 4PSK, 8PSK, 16QAM, 64QAM, 256QAM, 4096QAM.

FIG. 18 parameters 121 are used to calculate the MAP trellis and ML demodulation performance. MAP parameter values 122 are selected to provide best values for b while avoiding excessive numbers of layers n_(p) which increase the sidelobe induced demodulation loss and increase the required synchronization accuracy to support the n_(p). Plot point 5 lists two alternative candidates to achieve b=24 Bps/Hz which are: 64QAM requires n_(p)=6 layers and E_(b)/N_(o)=17.3 dB, and 256QAM requires a lower n_(p)=4 layers and a slightly higher E_(b)/N_(o)18.2 dB. The 256QAM appears to be the preferred candidate since the n_(p) is somewhat lower. ML parameter values 1,2,3 for the n_(s)=3 data symbol group assume n_(p)=4 QLM layers and the resulting performance in FIG. 16,17 is considerable lower than the MAP performance in FIG. 16,17. ML performance approaches the MAP performance in FIG. 16,17 for a n_(s)=4 data symbol group with increased number of QLM layers to n_(p)=6.

FIG. 19 is a proof-of-concept demonstration that QLM can provide a data rate which is substantially faster than the Shannon rate in (1), (11) with n_(p)=1, and in FIG. 17, 18. We use as a reference the performance for 256QAM at the Shannon rate which is C/W=6 Bps/Hz in FIG. 17 for E_(b)/N_(o)=10.21 dB which is the highest data rate mode for cellular communications and for almost all of the other communications links and applications. We observe that starting with 16QAM at the Shannon rate which is C/W=3 Bps/Hz for E_(b)/N_(o)=3.68 dB we can increase the data symbol rate to n_(p)=2, 3 times the 3 Bps/Hz whereupon the QLM performance is C/W=6, 9 Bps/Hz respectively for E_(b)/N_(o)=6.69, 8.95 dB which C/W is observed to be substantially above the Shannon bound or rate. Likewise it is observed that starting with 64QAM at the Shannon rate which is C/W=4 Bps/Hz for E_(b)/N_(o)=5.74 dB we can increase the data symbol rate to n_(p)=2, 3 times the 4 Bps/Hz whereupon the QLM performance is C/W=8, 12 Bps/Hz respectively for E_(b)/N_(o)=8.75, 11.01 dB which C/W is observed to be substantially above the Shannon bound or rate. This completes the proof-of-concept demonstration which complements the proof-of-concept in FIG. 17, 18.

FIG. 20 is an OFDM transmitter block diagram modified to support OFDM QLM with frequency offsets as the differentiating parameter to increase the symbol transmission rate from 1/T_(s) to the QLM rate n_(p)/T_(s) and with an increase in transmitter power to support this increased data rate. FIG. 20 differs from the IEEE WiFi 802.16 standard in the deletion of the 4.0 μs delay between the FFT output NT and the Tx FFT output NT_(s)=NT+4.0 μs, and in the use of the complete FFT band for transmission, in order to apply to all of the IEEE standard OFDM systems. Ideal OFDM modulates N input data symbols at the sample rate 1/T_(s) over the time interval NT_(s) with an N-point inverse fast fourier transform FFT⁻¹ to generated N harmonic waveforms e^j2πkn/N with each modulated by the corresponding data symbol wherein the normalized frequencies k=0,1, . . . N−1 correspond to channels 0,1, . . . , N−1, “j”=√(−1), “π=pi”, and “n” is a time index,. Data symbol output rates are 1/NT_(s) per channel and the N channels have a total symbol rate equal to N/NT_(s)=1/T_(s)=B=(Nyquist sample rate). Signal processing starts with the stream of user input data words (d_(k)} 46 with k indexed over the words. Frame processor 47 accepts these data words and performs turbo error correction encoding, error detection cyclic redundant encoding CRC, frame formatting, and passes the outputs to the symbol encoder 48 which encodes the frame data words into data symbols for handover to the OFDMA QLM signal processing. QLM transmits in parallel N Rx data symbols for each of the n_(p) FFT⁻¹ signal processing steams. Each set of Rx N data symbols are offset in frequency by 0, Δk, 2Δk, . . . , (n_(p)−1)Δk with Δk=1/n_(p) using the normalized frequency index k and are implemented in 49 by the frequency translation operator with FFT⁻¹ time sample index n. Following this frequency translation and FFT⁻¹ signal processing, the output streams of the OFDMA encoded symbols for the n_(p) frequency offsets are summed 51 and waveform encoded. The output stream of up-sampled complex baseband signal samples 52 {z(t_(i))} at the digital sample times t_(i) with digitization index i, is handed over to the digital-to-analog converter DAC, and the DAC output analog signal z(t) is single sideband SSB upconverted 52 to RF and transmitted as the analog signal v(t) wherein v(t) is the real part of the complex baseband signal z(t) at the RF frequency. Non-ideal OFDMA has a separation interval between contiguous FFT⁻¹ data blocks to allow for timing offsets and the rise and fall times of the channelization filter prior to the FFT⁻¹ processing.

FIG. 21 is a OFDM receiver block diagram modified to support OFDM QLM from the OFDM QLM transmitter in FIG. 20. Receive signal processing for QLM demodulation starts with the wavefronts 54 incident at the receiver antenna for the n_(u) users u=1, . . . , n_(u)≦N_(c) which are combined by addition in the antenna to form the receive Rx signal {circumflex over (v)}(t) at the antenna output 55 where {circumflex over (v)}(t) is an estimate of the Tx signal v(t) 52 in FIG. 20 that is received with errors in time Δt, frequency Δf, and phase Δθ. This Rx signal {circumflex over (v)}(t) is amplified and downconverted to baseband by the analog front end 56, synchronized (synch.) in time t and frequency f, waveform removed to detect the Rx QLM signal at the QLM symbol rate, inphase and quadrature detected (I/Q), and analog-to-digital ADC converted 57. ADC output signal is demultiplexed into n_(p) parallel signals 58 which are offset in frequency by 0, −Δk, −2Δk, . . . , −(n_(p)−1)Δk wherein Δk=1/n_(p) and processed by the FFT's. Outputs are trellis decoded 59 and further processed 60,61 to recover estimates of the Tx data d_(k) with k indexed over the data words.

FIG. 22 is an ideal LTE SC-OFDM transmitter block diagram modified to support SC-OFDM QLM with timing offsets as the differentiating parameter to increase the symbol transmission rate from 1/T_(s) to the QLM rate n_(p)/T_(s) and with an increase in transmitter power to support this increased data rate. FIG. 22 is ideal in the assumption of the use of orthogonal Wavelet multiplexing OWDM disclosed in patent U.S. Pat. No. 7,376,688 wherein the Wavelets provide ideal communications performance in that the symbol rates 1/T_(s) are equal to the bandwidth B=1/T_(s) and they are orthogonal in frequency offsets 1/T_(s) and time offsets T_(s), and the transmission is over a single channel and bandwidth in order to apply to all of the IEEE standard SC-OFDM LTE systems. Obviously, the transmitter block diagram applies to other waveforms equivalent to the Wavelet. Signal processing starts with the stream of user input data words (d_(k)} 111 with k indexed over the words. Frame processor 112 accepts these data words and performs the turbo error correction encoding, error detection cyclic redundant encoding CRC, frame formatting, and passes the outputs to the symbol encoder 113 which encodes the frame data words into data symbols for handover to the OWDM QLM transmit signal processing. The n_(p) time delays 0, ΔT_(s), 2ΔT_(s), 3ΔT_(s), . . . , (n_(p)−1)ΔT_(s) wherein ΔT_(s)=T_(s)/n_(p), are performed 114 and the output streams of the OWDM waveform encoded 115 symbols for the n_(p) time delays are summed 116 and passband waveform encoded and the up-sampled output stream of complex baseband signal samples 117 {z(t_(i))} at the digital sample times t_(i) with digitization index i, is handed over to the DAC and the DAC output analog signal z(t) is single sideband SSB upconverted 117 to RF and transmitted as the analog signal v(t) wherein v(t) is the real part of the complex baseband signal z(t) at the RF frequency. It is convenient to refer to the OFDM LTE communications as LTE(OFDM).

FIG. 23 is a LTE SC-OFDM QLM receiver block diagram modified to support SC-OFDM QLM from the transmitter in FIG. 22. Receive signal processing for QLM demodulation starts with the wavefronts 131 incident at the receiver antenna for the n_(u) users u=1, . . . , n_(u)≦N_(c) which are combined by addition in the antenna to form the receive Rx signal {circumflex over (v)}(t) at the antenna output 132 where {circumflex over (v)}(t) is an estimate of the Tx signal v(t) 117 in FIG. 22 that is received with errors in time Δt, frequency Δf, and phase Δθ. This Rx signal {circumflex over (v)}(t) is amplified and downconverted to baseband by the analog front end 133, synchronized (synch.) in time t and frequency f, waveform removed to detect the Rx QLM signal at the QLM symbol rate, inphase and quadrature detected (I/Q) and analog-to-digital ADC converted 134. ADC output signal is demultiplexed into n_(p) parallel signals 135 which are offset in time by 0, ΔT_(s), 2ΔT_(s), . . . , (n_(p)−1)ΔT_(s) and processed by the OWDMA decoders. Outputs are trellis decoded 136 with an algorithm comparable to the algorithm defined in FIG. 9. Outputs are further processed 137,138 to recover estimates of the Tx data d_(k) wherein k is indexed over the data words.

FIG. 24 illustrates a representative segment of current cellular communications networks with a schematic layout of part of a cellular network which depicts cells 1,2,3,4 that partition this portion of the area coverage of the network, depicts a user (network user) 5 located within a cell with forward and reverse communications links 6 with the cell-site base station (access point/hub) 7, depicts the base station communication links 8 with the MSC (mobile switching center) or the WSC (wireless switching center) 9, and depicts the MSC/WSC communication links with another base station (access point/hub) 17, with another MSC/WSC 16, and with external elements 10,11,12,13,14,15. One or more base stations are assigned to each cell or multiple cells or sectors of cells depending on the application. One of the base stations 9 in the network serves as the MSC/WSC which is the network system controller and switching and routing center that controls all of user timing, synchronization, and traffic in the network and with all external interfaces including other MSC's. External interfaces could include satellite 10, PSTN (public switched telephone network) 11, LAN (local area network) 12, PAN (personal area network) 13, UWB (ultra-wideband network) 14, and optical networks 15. As illustrated in the figure, base station 7 is the nominal cell-site station for cells i−2, i−1, i, i+1 identified as 1,2,3,4, which means it is intended to service these cells with overlapping coverage from other base stations (access points/hubs). The cell topology and coverage depicted in the figure are intended to be illustrative and the actual cells could be overlapping and of differing shapes. Cells can be sub-divided into sectors. Not shown are possible subdivision of the cells into sectors and/or combining the cells into sectors. Each user in a cell or sector communicates with a base station which should be the one with the strongest signal and with available capacity. When mobile users cross over to other cells and/or are near the cell boundary a soft handover scheme is employed for CDMA in which a new cell-site base station is assigned to the user while the old cell-site base station continues to service the user for as long as required by the signal strength.

FIG. 25 depicts a MIMO scenario for the cellular networks in FIG. 24 (2) using the multiple access OFDM QLM and SC-OFDM QLM communications encoding of the encoded data symbols transmitting signal information, followed by space-time coding H₀,C₀,L as disclosed in patent patent U.S. Pat. No. 7,680,211 for the communications architecture consisting of N Tx antennas 106 communicating through a scintillated, dispersive, fading, and multipath links 107 to M Rx antennas 108 for the forward communications links for the cellular communications 6 from the access point or hub 7 to the network user 5 in FIG. 24. Transmission starts with the input Tx data d 101 which is encoded, interleaved, formatted, and data symbol encoded 102 followed by the OFDM QLM and SC-OFDM QLM encoding 103 to generate the input signal (symbol) vector X 104 whose elements are the encoded data symbols from 103. This encoding includes the Tx signal processing required to support Rx signal processing equalization such as frequency domain equalization (FDE) and Wavelet domain equalization (WDE) to correct for scintillation, dispersion, fading, and multipath. Tx symbol vector X is space-time encoded 105 using H₀/H,C₀, L where (H₀/H) is equal to H₀ with the elements of H removed since they are generated by the Tx-to-Rx communication links, and handed over to the Tx antenna subsystem which performs the digital-to-analog conversion, symbol waveform modulation, single-side-band upconversion, power amplification, and transmission by the antenna elements for each of the N Tx transmissions by the access point or hub. Rx signals received by the network user antennas are amplified, filtered, downconverted, detection filtered and analog-to-digital converted to recover the Rx symbol vector Y 109. Space-time decoding 110 defined in scenarios 1,2,3,4 is implemented to recover the ML estimates {circumflex over (x)} 111 of X which are processed by the OFDM QLM and SC-OFDM QLM decoding and equalization 112 followed by deinterleaving and turbo decoding 113 to recover the ML estimates {circumflex over (d)} 114 of the Tx data d 101.

A multi-scale MS code can be implemented with modest complexity in order to improve the bit-error-rate BER performance of OFDM/LTE(OWDM), WiFi, WiMax, SC-OFDM QLM by spreading each Tx data symbol over the data band and over the data packet. This technique has been disclosed in U.S. Pat. No. 7,907,512. Jensen's inequality from mathematical statistics proves that this uniform spreading of the transmit signals using MS provides the best communications BER performance.

This patent covers the plurality of everything related to QLM generation for WiFi, WiMax, LTE, OFDM, and SC-OFDM waveforms, QLM demodulation for WiFi, WiMax, LTE, OFDM, and SC-OFDM waveforms, to QLM generation and data recovery of QLM, to the corresponding bounds on QLM, to all QLM inclusive of theory, teaching, examples, practice, and to implementations for related technologies. The representative trellis and ML algorithms for QLM demodulation are examples to illustrate the methodology and validate the performance and are representative of all QLM demodulation algorithms including all maximum likelihood ML architectures, maximum a posteriori MAP trellis and ML trellis architectures, maximum a priori, finite field techniques, direct and iterative estimation techniques, trellis symbol and iterative trellis symbol and with/without simplifications, trellis bit and iterative trellis bit and with/without simplifications and with/without bit error correction coding, sequential relaxation, equalization, super-resoluion, and all other related algorithms whose principal function is to recover estimates of the Tx symbols and information for QLM layered modulation and QLM increase in the data symbol rates as well as data recovery related to QLM and the QLM bounds.

This patent covers the plurality of everything related to QLM which has been disclosed in this and previous patents to consist of all methods and means to increase the communications data rates while introducing scaling or other possible performance improvements or combinations thereof to compensate for the impact of ISI and other communications link imperfections in reducing the available signal energy for demodulation and decoding to recover the communications information. The example applications in this invention disclosure have used a restricted version of QLM in order to derive the data rate performance to compare with the Shannon rate. 

What is claimed is:
 1. A method for implementation of a new Quadrature Layered Modulation (QLM) architecture for communications over the same frequency bandwidth of a carrier frequency, said method comprising the steps: introducing a new scaling parameter £ for the communications data symbol metrics E_(b)/N_(o), SNR to correct the n_(p) scaling for n_(p)>2 to compensate for the loss in the detectable signal energy caused by the interference of the overlapping signals, where said scaling parameter is a function of n_(p), the data symbol modulation, and the demodulation algorithm, and n_(p) is the increase in the data symbol rate using QLM and equivalently is the number of QLM channels when considering QLM to be a layered communications link, generating a scaled £n_(p)E_(b)/N_(o), £n_(p) ²SNR communications signal over the frequency bandwidth at the carrier frequency by modulating a stream of data symbols with a waveform at a n_(p)/T_(s) symbol rate wherein “T_(s)” is the time interval between contiguous symbols for the original channel prior to implementation of QLM, E_(b) is the energy per bit, N_(o) is the noise power density, and SNR is the signal to noise power ratio, transmitting and receiving said scaled communications signal over a QLM communications link consisting of a modulated stream of data symbols with a waveform at a n_(p)/T_(s) symbol rate, recovering data symbols of the scaled communications signals in a receiver using a demodulation algorithm, and combining said algorithm with error correction decoding to recover the transmitted information; whereby time offsets between contiguous data pulses is the differentiating parameter to enable the QLM data symbol rate which is independent of the Nyquist rate, to be demodulated.
 2. A method for implementation of a new Quadrature Layered Modulation (QLM) for communications over the same frequency bandwidth of a carrier frequency, said method comprising the steps: introducing a scaling parameter £ for the communications data symbol metrics E_(b)/N_(o), SNR to correct the n_(p) scaling for n_(p)>2 to compensate for the loss in the detectable signal energy caused by the interference of the overlapping signals, where said scaling parameter is a function of n_(p), the data symbol modulation, and the demodulation algorithm, and n_(p) is the increase in the data symbol rate using QLM and equivalently is the number of QLM channels when considering QLM to be a layered communications link, generating a first scaled £n_(p)E_(b)/N_(o), £n_(p) ²SNR communications signal over the frequency bandwidth at the carrier frequency by modulating a stream of data symbols with a waveform at a 1/T_(s) symbol rate wherein “T_(s)” is the time interval between contiguous symbols for the original channel prior to implementation of QLM, E_(b) is the energy per bit, N_(o) is the noise power density, and SNR is the signal to noise power ratio, generating a second scaled communications signal over the same frequency bandwidth at the same carrier frequency for a second channel by modulating a second stream of data symbols with the same waveform at the same symbol rate as the first stream of data symbols and with a time offset ΔT_(s) equal to ΔT_(s)=T_(s)/n_(p) wherein “n_(p)” is the number of QLM channels in said frequency bandwidth, for any additional channels, continuing generation of scaled communication signals over the same frequency bandwidth at the same carrier frequency by modulating additional streams of data symbols with the same waveform at the same data symbol rate as the first and second streams of data symbols, with time offsets increasing in each communication signal in increments of ΔT_(s)=T_(s)/n_(p) until the n_(p) signals are generated for n_(p) QLM channels, transmitting and receiving said scaled communications signals over a QLM communications link consisting of the n_(p) QLM channels, recovering data symbols of the scaled communications signals in a receiver using a demodulation algorithm, and combining said algorithm with error correction decoding to recover the transmitted information; whereby time offsets have been used as a differentiating parameter to enable QLM channels of communications over the same frequency bandwidth at the same carrier frequency with a data symbol rate independent of the Nyquist rate, to be demodulated.
 3. A method for implementation of a new Quadrature Layered Modulation (QLM) for communications over the same frequency band of a orthogonal frequency division multiplexed (OFDM) signal, said method comprising the steps: introducing a scaling parameter £ for the communications data symbol metrics E_(b)/N_(o), SNR to correct the n_(p) scaling for n_(p)>2 to compensate for the loss in the detectable signal energy caused by the interference of the overlapping signals, where said scaling parameter is a function of n_(p), the data symbol modulation, and the demodulation algorithm, and n_(p) is the increase in the data symbol rate using QLM and equivalently is the number of QLM channels when considering QLM to be a layered communications link, generating a first set of scaled £n_(p)E_(b)/N_(o), £n_(p) ²SNR communications signals at the available orthogonal frequency slots of the OFDM signal packet by modulating a first stream of data symbols at the orthogonal frequency slots spaced at 1/T_(s) intervals wherein “T_(s)” is the length of the of each QLM signal in an orthogonal frequency slot and is the length of the OFDM signal packet, E_(b) is the energy per bit, N_(o) is the noise power density, and SNR is the signal to noise power ratio, generating a second set of scaled communications signals at the frequency offset Δf=1/n_(p)T_(s) from the available orthogonal frequency slots of the OFDM signal packet by modulating said second stream of data symbols at these offset frequency slots at 1/T_(s) intervals, for any additional sets of channels, continuing generation of scaled communications signals with frequency offsets increasing in increments of ΔT_(s)=T_(s)/n_(p) from the available orthogonal frequency slots of the OFDM signal packet by modulating said stream of data symbols at these offset frequency slots at 1/T_(s) intervals until the n_(p) sets of signals are generated for n_(p) sets of QLM channels, transmitting and receiving said scaled communications signals over a QLM communications link consisting of the n_(p) sets of QLM channels, recovering data symbols of the scaled communications signals in a receiver using a demodulation algorithm, and combining said algorithm with error correction decoding to recover the transmitted information; whereby frequency offsets have been used as a differentiating parameter to enable the QLM parallel sets of channels of communications over the same frequency bandwidth at the same carrier frequency with a data symbol rate independent of the Nyquist rate, to be demodulated.
 4. The method of claim 1, 2, or 3 wherein the communication signals have the following properties: scaling parameter £ for the communications data symbol metrics E_(b)/N_(o), SNR corrects the n_(p) scaling for n_(p)>2 to compensate for the loss in the detectable signal energy caused by the interference of the overlapping signals, maximum normalized data rate (C/W) which is the communications efficiency in bps/Hz=b/s/Hz is defined by equation ${C/W} = {\max\limits_{np}\left\lbrack {n_{p}{\log_{2}\left( {1 + {{{SNR}/\pounds}\; n_{p}^{2}}} \right)}} \right\rbrack}$ wherein the maximum “max” is with respect to n_(p), “log₂” is the logarithm to the base 2, “W” is the frequency band in Hz, “SNR” is the ratio signal-to-noise over “W”, “bps” is the number of bits per second of the capacity “C” of the communications link, minimum signal-to-noise ratio per bit “E_(b)/N_(o)” is defined by equation ${\min\left\{ {E_{b}/N_{o}} \right\}} = {\min\limits_{np}\left\lbrack {\left( {{2\hat{}\left( {\left( {C/W} \right)/n_{p}} \right)} - 1} \right)\left( {\pounds\;{n_{p}^{2}/\left( {C/W} \right)}} \right)} \right\rbrack}$ wherein minimum “min” is the minimum with respect to n_(p), SNR=(E_(b)/N_(o))/(C/W), “E_(b)” is the energy per bit of the data symbol, “N_(o)” is the noise power density of the data symbol, QLM performance metric C/W is defined by equation C/W=n _(p) b wherein “b” is the number of information bits in the QLM data symbol, lower bound on the QLM performance metric E_(b)/N_(o) is defined by equation E _(b) /N _(o) =£n _(p)(2^−1)/b, and the lower bound on the QLM performance metric SNR is defined by equation SNR=£n _(p) ²(2^b−1) wherein these performance bounds apply to communications receiver demodulation performance of a QLM communications link consisting of n_(p) QLM channels or equivalently a QLM data symbol rate increase n_(p).
 5. The method of claim 1, 2, or 3 for implementation of QLM communications links, further comprising: introducing a scaling parameter £ for the communications data symbol metrics E_(b)/N_(o), SNR to correct the n_(p) scaling for n_(p)>2 to compensate for the loss in the detectable signal energy caused by the interference of the overlapping signals, where said scaling parameter is a function of n_(p), the data symbol modulation, and the demodulation algorithm, and n_(p) is the increase in the data symbol rate using QLM and equivalently is the number of QLM channels when considering QLM to be a layered communications link, constructing QLM channels of scaled £n_(p)E_(b)/N_(o), £n_(p) ²SNR communications E_(b)/N_(o), SNR for a QLM communications link over a frequency bandwidth of a carrier frequency with a data symbol rate independent of the Nyquist rate and E_(b) is the energy per bit, N_(o) is the noise power density, and SNR is the signal to noise power ratio, using one or more scaled data symbol modulations for the QLM communications channels, using one or more differentiating parameters to enable the channel signals of the QLM link to be demodulated, and implementing said demodulation using a trellis algorithm or maximum likelihood packet algorithm or a plurality of detection algorithms combined with error correction decoding algorithms.
 6. A Quadrature Layered Modulation (QLM) method for communications at faster-than-Nyquist (FTN) rate transmission of the data symbols, said method having the following properties: introducing a scaling parameter £ for the communications data symbol metrics E_(b)/N_(o), SNR to correct the n_(p) scaling for n_(p)>2 to compensate for the loss in the detectable signal energy caused by the interference of the overlapping signals, where said scaling parameter is a function of n_(p), said data symbol, said data symbol modulation, and the demodulation algorithm, and n_(p) is the increase in the data symbol rate using QLM and equivalently is the number of QLM channels when considering QLM to be a layered communications link, considering the data symbols at a FTN transmission rate to be layers of successive communications channels each at the Nyquist rate and with the last channel at a rate which is less than or equal to the FTN rate, scaling the communications performance metrics £n_(p)E_(b)/N_(o), £n_(p) ²SNR to enable the communications to maintain the same error rate performance for all FTN rates, where E_(b) is the data pulse signal energy per information bit, and N_(o) is the noise power density, SNR=S/N is the signal-to-noise power ratio or equivalently the signal-to-noise energy ratio of the data pulse, S=E_(b)b/T_(s) is the signal power, N=N_(o)/T_(s) is the noise power, N_(o) is the noise power density, SNR=b E_(b)/N_(o) where “b” is the number of information bits in the encoded data symbol, n_(p) is the increase in the data symbol rate using QLM, observing that QLM operates the link at a faster than Nyquist (FTN) rate since the Nyquist rate does not place any constraints on the information in a given bandwidth and therefore enables QLM to transmit more information than supported by the Shannon bound, demodulating the QLM communications at the data symbol rate to recover the transmitted data symbols, and combining said demodulation algorithm with error correction decoding to recover the transmitted information. 